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Abstract 

We show that the Borel resummed perturbative static potential at Nj = converges weU, and 
is in a remarkable agreement with the quenched lattice calculation at distances 1/r ^ 660 MeV. 
This shows that Borel resummation is very good at handling the renormalon in the static potential 
(and in the pole mass), and allows one to use the pole mass in perturbative calculation of heavy 
quark physics. 
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I. INTRODUCTION 



The asymptotic freedom of quantum chromodynamics (QCD) allows one to calculate 
the short distance physics accurately using perturbation. Unexpectedly, however, the per- 
turbative expansion of the static potential between a quark-antiquark pair does not show 
a convergence even at very short distances (see Fig. |1|). Moreover, no agreement is seen 
with the accurate lattice calculations of the static potential. This led to a suggestion of 
nonperturbative linear potential at short distance which, if proven true, would violate 
the expectation of the operator product expansion (OPE) that the nonperturbative effect 
at short distance is at most a quadratic potential. 

On the other hand, the bad convergence behavior of the perturbative expansion of the 
potential is well understood to be caused by the infrared (IR) renormalon which induces a 
constant nonperturbative effect proportional to Aqcd 0- This prompted several approaches 
to the problem. One is based on the observation that the force between a pair of static quarks 
is free from the leading renormalon. The potential obtained by integrating the force calcu- 
lated in perturbation indeed agrees quite well at short distance with the lattice potential 
P], up to an r independent constant. Another approach is the renormalon subtracted (RS) 
scheme [Q, in which one subtracts order by order the renormalon contribution from the 
perturbative potential. The potential calculated in this way also shows an improved conver- 
gence and agreement with the lattice potential. Another idea is to employ the cancellation 
of the renormalons in the static potential and the pole mass of the heavy quark |^. By 
expanding the pole mass and the static potential of a color singlet quarkonium in the run- 
ning coupling Osifi) and a short distance mass m(/i) one can avoid the renormalon problem, 
and indeed such an expansion shows an improved convergence [0) i|) 0- 

In this paper we show a more direct approach to the problem is possible via the Borel 
resummation of the perturbative potential. Since one might believe that the presence of 
an IR renormalon makes Borel resummation impossible, we state in advance that it is 
perfectly possible in this case. An IR renormalon in Borel resummation merely demands a 
corresponding nonperturbative effect, and since in this case it is a constant, the r dependence 
of the potential can be resummed with no difficulty. Moreover, this renormalon caused 
nonperturbative effect could be computed in the framework introduced in |TD[ where the 
nonperturbative effect is determined based on its conjectured analyticity in the complex 
coupling plane. An obvious advantage of the direct resummation is that the normalization of 
the potential can be fixed. In the approaches based on the renormalon cancellation/absence 
the potential can be fixed only up to an r independent constant. 

As we shall see the Borel resummed potential at short distance converges quickly, and 
agrees remarkably well with the lattice calculation, in fact better than any other approach 
introduced so far. The implication of this is significant. In the perturbative calculation of a 
heavy quark system one does not have to give up the pole mass in favor of a short distance 
mass to avoid the renormalon problem, and still can have a tight control on the perturbative 
expansion. 

Throughout the paper, unless stated otherwise, we consider pure QCD with no active 
quark flavors {Nf = 0), and the perturbative expansions considered are assumed to be in 
the MS scheme. As for the renormalon, we restrict our attention to the leading infrared 
renormalon that is closest to the origin in the Borel plane. 
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FIG. 1: The static potential at leading order (dotted), next-leading order (dashed), and next-next- 
leading order (solid). The data points denote lattice potential. 



II. BILOCAL EXPANSION OF THE BOREL TRANSFORM 

In general the perturbative expansion in weak coupling constant is an asymptotic expan- 
sion. When the large order behavior of the expansion is sign alternating like in 0^ theory 
it may be Borel resummed. However, when the expansion is of same sign at large orders 
Borel resummation demands a more careful treatment ITTI . In the case of the latter, one can 



first do Borel resummation at an unphysical negative coupling, at which the series is sign 
alternating, and then do analytic continuation in the complex coupling plane to the physical 
positive coupling. The Borel resummed amplitude obtained in such a way, however, turns 
out to have a cut along the positive real axis in the coupling plane, and consequently has an 
ambiguous imaginary part at a physical coupling. In Borel integration this imaginary part 
arises precisely from the infrared (IR) renormalon singularity of the Borel transform on the 
integration contour. This unphysical, ambiguous imaginary part then must be canceled by 
the nonperturbative effect corresponding to the renormalon. For further details we refer to 

0. 



Thus the static inter-quark potential V(r)^, which has an IR renormalon, can be written 
as the sum of the Borel integration with a contour on the upper (or lower) half plane and 
the nonperturbative effect [0, 

V\r, a,(l/r) ±ie] = — e-''/'^°"=(^/''V(6) db + V^plr, a,(l/r) ± ie] (1) 

rpo J0±ie 

where (3o is the one loop coefficient of the QCD (3 function, 

/3{as) = /i^^"* 



a^(/3o + /3ia, + /32a,' + ...), (2) 



^ Because of its infrared sensitivity the static potential is dependent on the ultrasoft factorization scale 
beginning at NNNLO |l^, however, to the order we are concerned (NNLO) this can be ignored. 
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and V{b) is the Borel transform that is given by 



with Vn defined in the perturbative expansion of the potential, 



(3) 



+1 



(4) 



l^NP denotes the renormalon caused nonperturbative effect. Since the imaginary parts in the 
first term in Eq. (|lD and in Vnp, respectively, cancel, the potential can be written as 



V[r,as{l/r)] = ^Re 



ooitie 



0±ie 



-fe//3oas(l/r) 



V{b) db 



+ Re{V^Np[r,«.(l/r)±ze]} . (5) 



Since I^np is an r independent constant proportional to Aqcd we can ignore it as far as the 
r dependence of the potential is concerned. However, a discussion on its determination will 
be given later on. 

The cancellation of the imaginary parts in the integral term and in Eq. (|ip determines 
the renormalon singularity in the Borel transform V{b). By comparing the functional form 
of 



Vnp oc A' 



MS 



l-^(/52/3o-A')//5o'«.(l/r) + 



(6) 



with the imaginary part of the Borel integration term in (|^), one can see V{b) must have 
the singularity 



cv 



1 + ci(l - 2b) + C2(l - 26)^ + . . . + Analytic part 



(7) 



;i - 2b)^+^ 

where the "Analytic part" denotes terms analytic around b = 1/2. The constants u and q, 
which depend only on the coefficients of the /3 function, were first determined in [H, and 



can be computed up to C2 from the known four loop /? function [15 



(3t + 4(31(3,(32 - 2Mlh + - 2f3!) - 2(3^^^ 
32u{u-l)Pl 



Pt - P0P2 



C2 



(8) 



The residue cy becomes the normalization constant of the large order behavior of the 
expansion (P, and its exact value is not known, but it can be determined perturbatively using 
the method developed in |jl6|, |l^]. Once cy is known, we can combine the two expansions of 
the Borel transform (^) and at 6 = and at 6 = 1/2, respectively, to obtain an improved 
description of the Borel transform in the region between the origin and the renormalon 
location at 6 = 1/2. There are in principle an infinite number of ways to interpolate the 
two expansions, but here we shall take a simple one which turns out to suffice our purpose 
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very well. We write the Borel transform as a two point expansion, which we call a bilocal 
expansion ^: 



Vib) 



lim Vn,m(&) 

N,M^cx) 



N 



Cv 



M 



l + ^c,(l-26)* 



i=l 



(9) 



By demanding that this bilocal expansion reproduce the expansion (^ around the origin 
the coefficients hn can be determined in terms of Vn and q. This gives, for example, the 
first three coefficients as 



ho = Vq- cv(1 + ci + C2) , 

hi = Vi- 2cy/3o[l - C2 + z/(l + Ci + C2)] , 

h2 = V2- icvPo[2 + z/(3 + ci - C2) + z/2(l + ci + C2)] . 



(10) 



For the bilocal expansion to work it is essential to have the residue cy calculated in a good 
accuracy, which is the subject of the next section. 



III. RENORMALON RESIDUE 



The residue can be determined in perturbation using the method developed in [[Iq , \L% . It 
was shown in 0, [l^, ^ that the residue in the case of the static potential can be calculated 
quite accurately. For completeness, we repeat the calculation here, and in the meantime 
obtain an improved estimate. 

To compute cy we ffist consider the function 

R{b) = {1- 2by+''V{b) . (11) 

Then, 

cv = Ri^) . (12) 
R{b) has a cut, but is bounded, at 6 = 1/2, and thus we can write cy as a convergent series, 

cv = j:rn(-) , (13) 



n=0 



where r„ are the coefficients of the power expansion of R{b) at the origin. The ffist three 
Vn can be calculated from the known Vn up to next-next-leading order (NNLO) [^, |2^, p3| , 
and this gives 

Cv ~ -1.33333 + 0.49943 - 0.33844 = -1.17234 . (14) 
The convergence is not that rapid but the series is oscillating. An important observation 



made in [19| is that the reliability of this estimate can be checked by the mutual cancellation 



of the renormalons in the static potential and the pole mass. 



This was first introduced in ||l8| in a slightly different context. 
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In perturbation theory the heavy quark pole mass mpoie can be expanded as 



n+l 



n=0 



(15) 



where mjj^ [= nT'jj^^mj^)] denotes the MS mass. As in the case of the static potential the 
Borel resummed pole mass can be written as 



1 + — 



ooitie 



g-V/3o".(mMs)mp„i^(6) db 



+mNp[a,(mj^) ±ie\ , 
where the Borel transform mpoie(6) has the perturbative expansion 

°° p f bV 



(16) 



(17) 



and mNp denotes the renormalon induced nonperturbative effect. The renormalon ambiguity 
in the pole mass proportional to A^ gives rise to a renormalon singularity that has exactly 
the same form as Eq. (^ of the static potential, 



rhpoicib) 



l + Ci(l-26) + C2(l-26)^ 



+ Analytic part , 



;i - 2by+'' 

Now the cancellation of the renormalons in 2mpoie and V{r) ^ leads to 

Cy + 2Cm = . 



(18) 



(19) 



We shall now compute the residue following the computation of cy. Using the known 
coefficients up to NNLO [|^, |2^, |2^ of the expansion (|l^) we have 

Cm ^ 0.42441 + 0.17473 + 0.02289 = 0.62203 (20) 

This time the convergence is quite good. With the two computed values we now have 

Cy + 2c„ 



Cy — 2c„ 



0.02968, 



(21) 



which shows a remarkable cancellation of the two residues. This gives an assurance on the 
accuracy of the calculated residues. 

We shall now compute c^ in a slightly different way. As has been shown in solvable 
models |jlO[, the knowledge on the renormalon locations in the Borel plane can be used 
in improving the residue calculation. Since we are interested in the power expansion of 
R{b) around the origin, we can obtain in principle a better convergence by expanding it in 
a new complex plane in which it is smoother around the origin |]18|. This can be done by 
pushing the renormalon singularities save the first one away from the origin with a conformal 
mapping. Let us consider the mapping []T^ P^ ] 



vTTb 



w 



26/3 



VT+b+ Jl-2b/3 



(22) 
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which maps the first renormalon at b = 1/2 to w = Wq, where 



Wo = ^, (23) 

and all other renormalons (at b = —n and b = 1/2 + n where n = 1, 2, 3, . . .) onto the unit 
circle. 

Expanding at the origin to O(w^) and evaluating it at w = Wq we have a new 

estimate of Cm 

^ 0.42441 + 0.16774 + 0.03451 = 0.62667 , (24) 

which is quite close to the previous one (^). This stability is reassuring that our computation 
is accurate. 

Now we shall quantitatively estimate the error in the computed residue (pij). We do this 
by computing Cm using an estimated NNNLO coefficient of the expansion (P^. We first 
estimate the unknown NNNLO coefficient following the method developed in First, 
expand R[b{w)] to 0{w^) with included. This gives 

R[b{w)] = 0.42441 + 0.83872W + 0.86284^^ + (-129.2687 + 3.43505 pg)^^ • (25) 

Note that the ps-independent constant term in the coefficient of is much larger than 
the coefficients of the lower orders. It turns out this is a generic feature of an asymptotic 
expansion with rapidly growing coefficients, and it can be used in estimating higher order 
unknown coefficients. From the pattern of the known lower order terms it appears quite 
reasonable to assume that the fourth coefficient is bounded by 

1129.2687 - 3.43505p3| < 2 . (26) 

This gives an estimate on ps 

P3 = 37.6322 ±0.58223. (27) 
With this result we can repeat the computation of Cm in w plane, now at NNNLO, to obtain 

= 0.62667 ±0.02553. (28) 

We thus conclude the error in the computed residue ( ^41) is about 4%. 

For the numerical analysis in Sec. we use the exact relation (|l^) and the pole mass 
residue ([2^ ) to compute cy. Since the convergence in the calculation of the pole mass residue 
is better than that of the potential, we would have a more accurate value this way. We thus 
have 

Cy = -1.25334 ±0.05106. (29) 



IV. DETERMINATION OF THE NONPERTURBATIVE EFFECT 

In this section we give an evaluation of the renormalon caused nonperturbative effect 
Vnp using the method developed in [1^. As mentioned in Sec. the role of Vnp in Borel 



resummation is to cancel the imaginary part arising from the renormalon singularity in the 
Borel integration of the static potential. This means that in principle the imaginary part 
of is calculable from perturbation theory. However, the real part, which is the physical 
quantity, cannot be directly calculated perturbatively without a further input. 
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The method for computing the real part rehes on the analyticity of Vnp in the complex as 
plane. As mentioned, the potential obtained by Borel resumming the asymptotic expansion 
has a cut along the positive real axis in the as plane, and from this cut the imaginary part 
of the perturbative term, the integral term in Eq. arises. To cancel this imaginary part 
it is thus plausible to demand that VNp(r, a^) also have a cut only along the positive real 
axis in the coupling plane. This then relates the real part to the perturbatively calculable 
imaginary part (we refer to [|TU| for details). For convenience, we shall call this method of 
determining the nonperturbative effect (along with the Borel integration of the perturbation 
series) 'analytic Borel resummation (ABR)'. Some nonperturbative effects in solvable models 
were shown to be calculable in ABR |jTo|. 

For ABR to work it is essential to have the functional form of the nonperturbative effect 
beforehand. In the case of the static potential it is provided by the renormalization group 
equation. Since Vnp in the MS scheme should be a constant proportional to Aj^, where 



^ME=-[(3oasil/r)]-''e 



u-l/2l3oOs{l/r) 



1 /•a^l/'") 



1 1 

+ 



dx 



(30) 



we can write, by demanding Vnp have a cut only along the positive real axis. 



VNp[r,a,(l/r)] = ^[-as{l/T)re-"'^^^^^"^^ 



X exp 




-A 

(3{x) ' (3oal (3lx 



1 1 

+ 



dx 



(31) 



with C an undetermined real constant. Note that a cut can arise only from the prefactor in 
Eq. ( pOD with a noninteger v. Now the cancellation of the imaginary part in Vnp[?", as(l/r) ± 
ie\ with the corresponding imaginary part in the Borel integration term in Eq. (Q) fixes the 
constant C: 

C=^At^. (32) 



(2/5o)i+^ 



The real part of Vnp is then given by 



Re [VNp(a, ± le)] = ^^1+^^^^ cos(z/7r)Aj^ . 

With the calculated residue cy in Eq. (p9|), we find at A"/ = 

Re [VNp(a, ± ze)] = 0.477Aj^. 



(33) 



(34) 



In the numerical analysis in the next section we will combine this result with the Borel 
integration of the perturbative expansions. 

V. COMPARISON WITH LATTICE CALCULATION 



The static potential in lattice calculation is extracted from the Wilson line of a static 
quark- ant iquark pair, computed in Monte Carlo simulation. The recent calculations [^, |2^, 
|30| , |3l| employing large lattices up to 64^ achieved a remarkable accuracy, and can probe a 
short distance where perturbative QCD should be applicable. It is thus an ideal place where 
perturbative QCD can be compared with lattice calculations. 
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FIG. 2: The lattice potential vs the Borel resummed potential using Vo,2 (dotted), Vi^2 (dashed), 
and V2,2 (solid). 



As we mentioned in Introduction, the truncated power series of the perturbative expansion 
fails even at very short distance. We shall now see this problem can be cured by Borel 
resummation. 

The numerical integration of the Borel integral in Eq. (|^) can be done easily in w plane 
defined by the mapping (p^ . Using the Cauchy's theorem, the integration contour, for 
example, on the upper half plane in w plane can be deformed to a ray off the origin to the 
unit circle in the first quadrant. This trick allows to avoid the renormalon singularity on 
the integration contour, and makes the computation easy. For details we refer to |1T8| . 

For comparison with lattice calculation we take the accurate data of the recent computa- 
tion employing large lattices [2S]. All the dimensional quantities are in units of the Sommer 
scale ro (~ 0.5 fm) where tq in terms of Aj^ (^ 238 Me V) is determined in lattice 
computation to be 

roAy^ = 0.602(48) . (35) 

On the side of the perturbative potential, the Borel integration in Eq. (|) was done using 
the Borel transform Vo,2,^i,2, and V2,2 in the bilocal expansion (^. The coupling constant 
as{l/r) was computed by numerically solving Eq. (|30D employing the four loop /? function 
|T5| . Because of the divergent quark self energy the lattice potential is determined only up 
to an r independent constant, so we subtracted such a constant from the lattice data so that 
the lattice potential and the NNLO perturbative potential agree exactly at r/ro = 0.30798. 

The result is in Fig. |^. Notice the rapid convergence of the resummed potential at 
distances r ^ 0.6ro [~ (660MeV)~^], and the excellent agreement of the NNLO potential 
with the lattice data. The potential at leading order already fits the lattice values quite 
well. It is remarkable that perturbative QCD is applicable at distances as large as r = 
(660MeV)-^ 
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VI. DISCUSSION AND SUMMARY 



The first thing we can learn form our result is that in the static potential the leading 
renormalon is overwhelmingly dominant at short distances and there cannot be any signifi- 
cant nonperturbative effect other than that caused by the renormalon. As already observed 
in 0, ^, large linear potentials at short distances like those proposed in 0, |35[ are 
excluded. 

The rapid convergence of the perturbative potential in ABR allows one to use the pole 
mass in perturbative calculation of heavy quarkonium physics. Because of the bad conver- 
gence of the truncated power series of the static potential, there was a limit in the precision 
achievable with perturbative QCD in quarkonium physics PB| , P7| . But, it was soon realized 
that the cancellation of renormalons in the pole mass and the static potential can be used 
to alleviate the problem |^. Instead of using the pole mass directly, one can achieve an 
improved convergence by simultaneously expanding the pole mass and static potential in the 
heavy quark Hamiltonian in terms of the running coupling as{n) and a short distance mass 
like the MS mass Q . Although this approach avoids the renormalon problem, there could 
be large logs in the perturbative expansion which could in principle spoil the convergence. 
Since the expansion involves two far-separated scales, the heavy quark mass and 1/r mv, 
where v is the heavy quark velocity) large logs like ln(r/i) and/or ln(m//i) could survive for 
any choice of fi, which in practice is typically taken values in-between the two scales. With 
our resummation of the static potential, the convergence problem at short distance is solved, 
so the pole mass needs not be abandoned in favor of a short distance mass. Once the pole 
mass is extracted by comparing, say, a calculated quarkonium spectrum to an experimental 
value, the MS mass can be obtained from the pole mass by resumming the quark mass 
expansion (p!5| ) in ABR. Since the renormalon in the pole mass is essentially same as that 
in the static potential, we can expect a rapid convergence of the Borel resummation of the 
mass expansion, and we have checked that this is indeed the case. As an example, for the 
bottom quark {Nf = 4) with as(m^) = 0.22 the 'Borel resummed (BR)' mass mBR, which 
is defined as the real part of the integral term in Eq. (0), converges as 

mBR = m^{l + 0.15769 + 0.00409 - 0.00028) . (36) 

Notice the rapid convergence. The renormalon caused nonperturbative effect m^p in Eq. 
(|16|) can be determined in ABR, and its real part equals to — Re[VNp]/2 that is given in Eq. 
(^). An obvious advantage of the direct resummation of the renormalons is the separation of 
scales; The perturbative expansions for the pole mass and the static potential are resummed 
at their optimal scales /i = Tny^ and fi = 1/r, respectively, and there is no mixing of these 
scales as in the above implementation of renormalon cancellation using a short distance 
mass. The absence of large logs and the excellent convergence of the resummed mass and 
potential are expected to provide a new level of precision calculation for heavy quarkonium. 

It is worthwhile to mention that the nonperturbative effects V^p and rriNp may actually 
decouple completely from the quarkonium system. The renormalon cancellation between 
the pole mass and the static potential means that the ambiguous imaginary parts in these 
quantities cancel without the introduction of the nonperturbative effects. This implies that 
the nonperturbative effects are actually spurious, appearing only at an intermediate step in 
Borel resummation, and physical observables are completely independent of them. Specifi- 
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FIG. 3: The strong couplings obtained by employing the four loop (3 function (solid) and its [2/3] 
Fade approximant (dashed). 

cally, we may write the Hamiltonian of a heavy quarkonium system as 

— >2 

H = 2mpoie + ^ + V[r, a.(l/r)] . (37) 

^pole 

Putting 

"^poie = "^BR[mj^, as(mj^)] + Re[mNp] , 
V[r, «.(l/r)] = VbkV, Ois{l/r)] + Re[l^Np] , (38) 

where the BR potential Vbr denotes the real part of the integral term in Eq. (|^), and using 
the cancellation of 2Re[mNp] with Re[VNp] in ABR,^ we can write H in terms of the BR 
quantities only: 

H = 2mBR + ^ + VBR[r, «s(l/r)] + O (p^Re[m^p] / ml^) . (39) 

The remaining dependence on the nonperturbative effect suppressed by an inverse power of 
the quark mass is expected to cancel when higher order terms in quark mass expansion of 
the Hamiltonian are taken into account. This shows that the Hamiltonian in BR scheme 
is formally same as that in the on-shell scheme with the on-shell quantities mpoie and V{r) 
replaced by the corresponding BR quantities. Thus for physical observables the specific form 
of the nonperturbative effects are not necessary. 

The perturbative potential and the lattice values in Fig. |^ begin to deviate at r ^ 
0.6ro, which we regard as the failure of the perturbative potential at these distances. It is 
interesting to observe that this deviation occurs approximately at the same position where 
the four loop f3 function fails. The couplings ^^(l/r) obtained by running with the four 
loop f3 function and its [2/3] Fade approximant, which differs from the former only at 



^ This cancellation is not automatic but a feature of ABR. 
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orders higher than four loop, are plotted in Fig. ^. Notice that they begin to deviate 
approximately at the same distance where the perturbative potential begins to fail. At 
r = 0.6ro [as{l/0.6ro) = 0.417] the P function has the expansion 

P = -0.152(1 + 0.308 + 0.143 + 0.097 + . . .) (40) 

which shows the convergence is quite slow at this distance. It seems the coupling grows too 
fast at these distances, since a more slowly growing coupling would fit the lattice data. This 
simultaneous deviations could be a coincidence, but a more plausible explanation would 
be that the failure of the (3 function at these distances results in an unreliable coupling, 
which then causes the deviation. The (3 function would not be all that fails the perturbative 
potential. Since there is a renormalon singularity at b = 3/2 the bilocal expansion (^ 
at a finite order would certainly fail around b ^ 3/2. This does not cause any serious 
problem at small couplings, but as the coupling increases this becomes problematic because 
the Borel integral in Eq. (^ receives a sizable contribution from the region far from the 
origin. By varying the upper bound of the integration in Eq. (^) one can easily check that 
the resummed potential at r ^ 0.6ro is indeed sensitive on the Borel transform at 6 ^ 3/2. 
This argument suggests that the applicability of the Borel resummed perturbative potential 
could be extended to larger distances once we have a better control over the (3 function and 
the Borel transform at such distances. 

Lastly, we note that the convergence problem of the truncated power series in the per- 
turbative potential is only one example, although a very conspicuous one, of the problem of 
the QCD expansions in general, especially, at low energies of a few GeVs. The problem was 
not so visible in these expansions, since many were considered at a fixed scale, not like the 
perturbative potential considered here where a continuum of scale is involved. Convention- 
ally, in the OPE approach, in these low energy expansions the physical quantity is organized 
as the sum of a truncated power series and power corrections. Any difference between the 
truncated power series and the (unknown) true value is swept over to the power correc- 
tions. Clearly, this approach fails in the static potential because the potential of the OPE 
approach is just the truncated power series plus an r independent constant, which we know 
has a bad convergence and disagrees with the lattice calculation. As already discussed more 



extensively in the Gross-Llewellyn Smith sum rule the solution to the problem is the 
Borel resummation that properly accounts for the renormalon. Without Borel resummation 
the bad convergence in the truncated power series results in wide fluctuations in the power 



corrections as the order of perturbation varies, which is observed in many cases. See ^ 
for some examples. 

To summarize, we have shown that the Borel resummation with a proper account of 
the renormalon singularity in the Borel plane can resolve the convergence problem of the 
perturbative static potential and the pole mass, and the potential obtained in such a way 
is in an excellent agreement with the lattice calculation. Consequently, any significant 
nonperturbative effect at short distance other than the renormalon effect is excluded, and 
the pole mass can be used in perturbative calculation of heavy quarkonium physics. The 
advantages of the direct resummation of the renormalons include rapid convergence of the 
summations and absence of large logs, and these can open a new level of precision calculation 
for heavy quarkonium. We also calculated in the framework of ABR the renormalon caused 
nonperturbative effects in the static potential and the pole mass. The resummation method 
developed here may be applied to the computation of heavy quarkonium spectra in an 
approach similar to that employed in [M, where the perturbative potential at short distance 
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is combined with the phenomcnological potential at large distance. Also it may be employed 
in the top threshold production. 
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